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A  multi-objective  and  multi-parameter  optimization  is  implemented  to  design  the  optimal  structure  of 
bismuth-telluride-based  TEG  (thermoelectric  generator)  module.  A  multi-physics  TEG  model  combining 
the  SCG  (simplified  conjugate-gradient)  algorithm  is  used  as  the  optimization  tool.  The  semiconductor 
pair  number,  leg  length,  and  base  area  ratio  of  semiconductor  columns  to  TEG  module  significantly  affect 
the  TEG  performance,  and  hence  are  all  incorporated  into  the  present  optimization  study.  A  single¬ 
objective  optimization  is  first  implemented  to  provide  input  parameters  for  the  multi-objective  opti¬ 
mization.  The  results  show  that  when  taking  the  output  power  as  the  single-objective  function,  the 
output  power  can  be  elevated  significantly  by  optimization  of  the  three  geometric  parameters  but  which 
also  accompanies  the  significant  reduction  in  the  conversion  efficiency.  The  same  result  also  occurs  when 
taking  the  conversion  efficiency  as  the  single-objective  function.  By  combining  the  output  power  and 
conversion  efficiency  with  a  weight  factor  as  the  multi-objective  function,  the  optimization  is  again 
implemented.  The  optimal  design  obtained  by  multi-objective  optimization  makes  a  proper  balance 
between  the  output  power  and  conversion  efficiency,  so  that  the  both  are  improved  simultaneously. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

TEGs  (thermoelectric  generators)  are  based  on  the  Seebeck  ef¬ 
fect  of  semiconductor  materials  to  convert  thermal  energy  to 
electricity  directly  [1,2].  In  recent  years,  thermoelectric  technology 
is  attracting  more  and  more  attention,  mainly  because  of  the 
following  two  aspects:  on  the  one  hand,  there  are  not  working 
fluids  or  other  moving  parts,  so  TEs  (thermoelectric  devices)  have 
many  good  features,  such  as  reliable  operation,  layout  flexibility, 
adaptability  and  other  characteristics  [3,4];  on  the  other  hand,  TEs 
do  not  produce  secondary  pollution  gases  such  as  carbon  dioxide  or 
other  unfriend  polluting  gases  in  the  progress  of  using  the  life  or 
industrial  waste  heat  for  electricity  generation  [5,6]. 

Previous  studies  show  that  the  TEG  performance,  including  the 
output  power  (P)  and  conversion  efficiency  (r;),  is  strongly  depen¬ 
dent  on  semiconductor  material  properties.  In  general,  the 
dimensionless  figure  of  merit  is  adopted  to  represent  the  influence 
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of  material  properties  on  the  TE  (thermoelectric  device)  perfor¬ 
mance,  which  is  defined  as  ZT  =  a2aTIXt  where  a  is  the  Seebeck 
coefficient,  a  is  the  electric  conductivity,  A  is  the  thermal  conduc¬ 
tivity,  and  T  is  the  absolute  temperature  at  which  the  properties  are 
measured  [7].  A  larger  ZT  represents  better  performance  of  TEs, 
therefore,  researchers  who  are  engaged  in  thermoelectric  materials 
are  committed  to  enhancing  ZT,  which  means  increasing  the  See¬ 
beck  coefficient  and  electric  conductivity,  and  at  the  same  time 
keeping  thermal  conductivity  values  as  low  as  possible  [8,9].  Min- 
nich  et  al.’s  research  [8]  showed  that  bulk  nanostructured  ther¬ 
moelectric  materials  have  the  most  promise  for  commercial  use 
due  to  their  plasticity  and  compatibility.  Tang  and  Zhang  [9]  found 
that  the  ZT  of  thermoelectric  material  FexCo4_xSbi2  (x  =  0-3)  can 
be  improved  significantly  by  filling  Ce,  Ba,  and  Y. 

Except  for  thermoelectric  materials,  the  impact  of  the  geometry 
on  the  TEs  performance  cannot  also  be  ignored.  TEs  geometry 
optimization  can  be  divided  into  two  categories:  one  is  the  single¬ 
parameter  optimization  [10-15],  and  the  other  is  multi-parameter 
optimization  using  various  inverse  problem  approaches  [16-19]. 
Up  to  now,  the  single-parameter  optimization  is  the  most  widely 
used  method.  For  example,  Jang  et  al.  [10  adopted  the  single¬ 
parameter  method  to  study  how  the  thickness  of  ceramic  plate, 
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Nomenclature 

Greek  letters 

a 

Seebeck  coefficient  (V  I<-1) 

P 

output  power  (W) 

a 

electric  conductivity  (S  m-1) 

V 

conversion  efficiency 

X 

thermal  conductivity  (W  m_1  K-1) 

COP 

coefficient  of  performance 

P 

electric  resistivity  (Q  m-1) 

ZT 

figure  of  merit 

A 

difference 

T 

temperature  (I<) 

P 

search  step  size 

f 

weight  factor 

TC 

search  direction 

C 

search  variable 

T 

base  area  ratio  of  semiconductor  columns  to  TEG 

^total 

area  of  TEG  module  base 

module 

J 

objective  function 

N 

semiconductor  pair  number 

Subscripts 

Tlpair 

base  area  of  TEG  unit 

opt 

optimal 

Tlpn 

cross-sectional  area  of  p-  or  n-type  semiconductor 

pn 

p-  or  n-type  semiconductor 

H 

height  (mm) 

Cu 

copper 

W 

width  (mm) 

k 

number  of  variables 

Abbreviation 

Superscript 

TEG 

thermoelectric  generator 

n 

number  of  iterations 

TEC 

thermoelectric  cooler 

foot  length  and  cross-sectional  area  of  the  p-  or  n-type  semi¬ 
conductor  affect  the  performance  of  a  micro  TEG  (thermoelectric 
generator)  unit.  They  found  that  for  each  of  the  three  geometry 
parameters,  there  is  always  a  specific  size  corresponding  to  the 
optimal  P  or  77,  however,  the  optimal  P  and  77  cannot  be  reached 
simultaneously.  Yu  et  al.  [11  ]  lay  emphasis  on  the  influences  of  the 
junction  temperature  difference  (AT2)  in  the  second  stage,  the 
length  of  thermocouples  (Li)  and  the  NUM  (number  of  thermo¬ 
couples)  in  the  first  stage  on  the  cooling  performance  of  a  two- 
stage  TEC  module.  They  included  that:  1 )  For  a  fixed  AT2,  the  de¬ 
creases  of  l\  and  NUM  can  improve  the  coefficient  of  perfor¬ 
mance  ;  2)  When  keeping  a  constant  coefficient  of  performance, 
the  decreases  of  L\  and  NUM  can  also  decrease  AT2,  and  thus  in¬ 
crease  the  cooling  capacity.  Chen  et  al.  [12]  proposed  a  new  cycle 
model  consisting  of  a  multi-couple  thermoelectric  device  and 
involving  several  key  irreversibilities  of  real  TEGs  to  optimize  the 
value  of  (Sp//p)n,  where  Sp  and  /p  are  the  cross-sectional  area  and  leg 
length  of  p-type  semiconductor,  n  is  the  number  of  couples.  For  the 
fixed  total  number  of  thermoelectric  elements  of  two-stage  ther¬ 
moelectric  refrigerator  driven  by  two-stage  thermoelectric  gener¬ 
ator,  Meng  et  al.  [13]  optimized  the  allocations  of  the 
thermoelectric  element  pairs  for  the  maximum  cooling  load  and 
coefficient  of  performance.  Yilbas  and  Sahin  [14]  examined  the 
influence  of  the  slenderness  ratio  and  the  external  load  parameter 
on  the  maximum  efficiency  of  a  TEG,  and  found  that  the  efficiency 
attains  high  values  for  the  slenderness  ratio  less  than  1  for  almost 
all  the  external  load  parameters  considered,  which  is  more  pro¬ 
nounced  for  the  large  values  of  the  external  load  parameter.  Opti¬ 
mizations  were  carried  out  for  the  two-stage  TE  coolers  with  two 
design  configurations,  i.e.,  the  pyramid-styled  and  cuboid-styled 
ones  by  Xuan  et  al.  [15],  the  optimum  number  ratio  of  TE  mod¬ 
ules  between  two  stages  for  the  first  design  and  the  optimum 
electric  current  ratio  between  two  stages  for  the  second  design 
were  finally  obtained. 

In  all  the  single-parameter  optimization  studies  10-15],  the 
optimal  value  for  a  specific  parameter  is  obtained  by  searching  this 
parameter  to  reach  the  optimal  TE  performance  while  keeping  the 
other  parameters  unchanged.  However,  because  all  the  parameters 
have  coupled  effects  on  the  TE  performance,  the  optimal  TE  per¬ 
formance  cannot  be  reached  only  by  searching  a  specific  parameter. 
Fortunately,  a  few  studies  [16-19]  have  made  attempt  to  develop 


TEC  (thermoelectric  cooler)  multi-parameter  optimization 
methods  by  combining  the  various  optimization  algorithms  and  the 
TEC  analytical/numerical  model.  Cheng  et  al.  combined  a  TEC 
model  and  a  genetic  algorithm  to  optimize  of  the  geometry  of 
single-  16]  and  double-stage  [17]  TEC.  In  their  study,  the  maximum 
cooling  capacity  was  defined  as  the  objective  function  under  the 
requirement  of  minimum  coefficient  of  performance  and  the 
constraint  of  maximum  material  cost.  The  leg  length,  the  leg  area 
and  the  number  of  legs  were  taken  as  search  variables  and  were 
optimized  simultaneously.  Rao  and  Patel  [18]  developed  a  modified 
version  of  the  teaching-learning-based  optimization  algorithm 
and  incorporated  this  algorithm  into  the  TEC  model  to  optimize  the 
performance  of  double-stage  TEC.  The  cold  stage  current,  hot  stage 
current,  and  number  ratio  of  thermoelectric  elements  between  the 
two  stages  were  three  search  variables.  Thermal  resistance  model 
was  adopted  as  the  direct  problem  solver  in  Refs.  [16-18].  Due  to 
grossly  simplifying  assumptions,  such  as  that  only  the  energy  bal¬ 
ance  equations  at  the  hot  and  cold  ends  were  solved,  only  constant 
or  temperature-averaged  properties  can  be  used,  the  current  den¬ 
sity  vector  was  simplified  as  a  constant  scalar,  etc.,  the  thermal 
resistance  model,  however,  cannot  predict  the  TE  performance 
accurately,  and  therefore  is  not  appropriate  as  the  direct  problem 
solver  for  TE  optimization.  Recently,  we  combined  a  multi-physics 
TEC  model  [19]  and  a  SCG  algorithm  to  optimize  the  cooling  ca¬ 
pacity  of  a  TEC  module.  A  value  of  COP  (coefficient  of  performance) 
that  was  not  lower  than  70%  COP  of  the  initial  TEC  geometry  was 
chosen  as  the  constraint  condition,  and  the  semiconductor  pair 
number,  the  leg  length,  and  the  base  area  ratio  of  semiconductor 
columns  to  TEC  were  taken  as  search  variables.  The  optimal  ge¬ 
ometry  corresponding  to  the  maximum  cooling  capacity  was  suc¬ 
cessfully  obtained,  however,  the  COP  of  the  optimal  geometry  was 
reduced  as  compared  with  the  initial  geometry. 

The  motivation  of  the  present  work  is  based  mainly  on  the 
following  three  points:  1)  the  effect  of  the  TEG  geometry  on  its 
performance  is  also  very  significant,  however,  the  TEG  multi¬ 
parameter  optimization  has  not  been  reported  so  far;  2)  as  two 
important  evaluation  indexes  for  the  TEG  performance,  P  and  77 
cannot  reach  the  extreme  value  simultaneously  when  using  the 
single-objective  optimization,  hence,  multi-objective  optimization 
needs  to  be  carried  out;  3)  it  is  well  known  that  the  direct  problem 
solver  is  very  critical  for  the  optimization,  due  to  the  limitation  of 
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the  thermal  resistance  model,  the  optimization  needs  to  adopt  the 
multi-physics  TEG  model  as  the  direct  problem  solver.  In  the  pre¬ 
sent  work,  the  TEG  multi-objective  optimization  approach  is  first 
developed,  which  is  composed  of  a  SCG  algorithm  and  a  multi¬ 
physics  TEG  model.  The  approach  is  then  used  to  optimize  a 
micro-TEG  module.  By  taking  the  weighted  average  of  P  and  p  as  the 
multi-objective  function,  three  search  variables,  semiconductor 
pair  number,  leg  length,  and  base  area  ratio  of  semiconductor  col¬ 
umns  to  TEG  module,  are  optimized  simultaneously  to  look  for  the 
optimal  TEG  performance. 

2.  Numerical  methods 

2.1.  Optimization  approach 

Typically,  an  inverse  problem  optimization  approach  includes  a 
direct  problem  solver  and  an  inverse  problem  solver.  To  ensure  the 
successful  optimization,  the  multi-physics  TEG  model  [20-22] 
developed  previously  by  our  group  is  adopted  as  the  direct  prob¬ 
lem  solver.  The  multi-physics  TEG  model  solves  the  temperature 
and  electric  potential  equations,  and  accounts  for  all  effects 
occurred  in  TEG,  including  Seebeck  effect,  Peltier  effect,  Thomson 
effect,  Joule  heating,  Fourier  heat  conduction,  and  temperature- 
dependent  properties.  The  TEG  model  was  proven  to  be  complete 
self-consistency  [20],  which  also  has  been  well  validated  by 
comparing  the  prediction  with  the  experimental  data  and  the 
analytical  solution  [21,22],  and  its  accuracy  is  greatly  improved  as 
compared  with  the  thermal  resistance  model  [23].  More  details  of 
the  model  can  be  found  in  Refs.  [20-22]. 

There  are  many  optimization  algorithms,  such  as  the  genetic 
algorithm,  conjugate-gradient  algorithm,  Newton-Raphson  algo¬ 
rithm,  adaptive  weighting  input  estimation  algorithm,  and  simu¬ 
lated  annealing  algorithm.  The  genetic  algorithm  and  simulated 
annealing  algorithm  search  in  the  direction  of  a  better  design,  but 
randomly  mutate  from  this  direction  in  case  there  might  be  a  better 
solution  elsewhere.  The  advantage  of  this  type  of  algorithms  is  that 
they  are  able  to  find  the  global  solution  among  multiple  local 
minima.  The  disadvantage  is  that  they  require  a  very  large  number 
of  design  evaluations  and  cannot  present  how  the  optimized  pa¬ 
rameters  approach  their  optimal  values.  Another  class  of  algo¬ 
rithms  is  known  as  gradient  algorithms.  These  algorithms  use 
information  about  how  a  change  in  each  design  variable  affects  the 
objective  function.  The  gradient  algorithms  such  as  conjugated- 
gradient  algorithms  calculate  the  partial  derivative  of  the  objec¬ 
tive  function  with  respect  to  each  design  variable  and  use  the  in¬ 
formation  to  select  design  variable  values  for  an  improved  design. 
The  advantage  of  gradient  algorithms  is  that  once  the  gradients  are 
calculated  they  are  efficient.  Therefore,  the  present  work  adopted 
the  conjugated-gradient  algorithm  as  the  optimization  algorithm. 
Determining  the  optimal  search  step  size  for  conventional 
conjugated-gradient  procedures  is  complex  and  time-consuming. 
Cheng  and  Chang  [24]  developed  a  simplified  version  of  the 
conjugated-gradient  algorithm  that  uses  a  constant  step  size  in 
each  iteration  and  applies  a  direct  numerical  differentiation 
method  to  determine  search  direction.  The  simplified  conjugated- 
gradient  algorithm  has  been  applied  successfully  to  optimize  the 
geometric  structures  of  proton  exchange  membrane  fuel  cells  [25], 
micro-channel  heat  sinks  [26],  and  TEC  19]  in  our  previous  studies. 
Therefore,  this  algorithm  is  chosen  as  the  inverse  problem  solver. 


unit  (Fig.  lb)  has  the  base  area  of  Apair  =  Atota\lN  and  consists  of 
three  metal  connectors,  a  p-type  semiconductor  and  a  n-type 
semiconductor.  The  p-  and  n-type  semiconductors  have  the  same 
dimension  with  the  height  of  Hpn  and  the  cross-sectional  area  of 
Apn  =  Wcu  x  Wcu.  The  height  of  the  copper  connector  is  Hcu  and  the 
width  is  Wcu.  The  metal  connector  is  made  of  copper,  the  p-  and  n- 
type  semiconductors  are  bismuth-telluride  materials.  Table  1  lists 
all  the  material  properties.  It  is  noted  that  the  Seebeck  coefficient, 
electrical  conductivity  and  thermal  conductivity  of  semiconductors 
are  all  temperature-dependent. 

For  the  TEG  module  with  a  fixed  Atotai.  once  the  semiconductor 
pair  number  (N),  the  leg  length  (Hpn),  and  the  base  area  ratio  of 
semiconductor  columns  to  TEG  module  (y  =  2NApn/Atotai  =  2Apn/ 
Apair)  are  specified,  the  TEG  module  geometry  will  be  determined 
uniquely.  The  effect  of  single  parameter  N,  Hpn  or  y  on  the  TEG 
module  performance  is  shown  in  Fig.  2.  Here  the  cold  end  tem¬ 
perature  is  assumed  to  be  300  K  with  the  temperature  difference  of 
AT  =  150  K.  N  =  40,  Hpn  =  0.5  mm  and  y  =  0.80  are  defined  as  the 
base-line  geometry.  In  the  single-parameter  analysis,  only  one 
parameter  of  interest  is  changed  while  the  other  two  parameters 
keep  the  values  of  the  base-line  geometry.  Fig.  2  shows  that  N,  Hpn 
and  y  all  significantly  affect  the  TEG  module  performance.  Conse¬ 
quently,  the  three  parameters  are  chosen  as  the  search  variables. 

Fig.  2  also  shows  that  for  the  output  power  the  optimal 
parameter  is  N  =  60,  Hpn  =  0.06  mm,  and  y  approaching  to  1, 
however,  it  becomes  N  =  55,  Hpn  =  0.68  mm,  and  y  =  0.60  for  the 
conversion  efficiency.  It  can  be  expected  that  when  the  output 


2.2.  TEG  module  geometry  and  optimized  parameters  Lc 

Fig.  la  is  the  schematic  of  a  TEG  module,  which  comprises  two 
bases  with  the  area  of  Atotai  =  4x4  mm2,  N  pairs  of  TEG  units  1  j 

connected  electrically  in  series  and  thermally  in  parallel.  Each  TEG  Fig.  l.  Schematics  of  a  teg  module  and  a  teg  unit:  (a)  teg  module;  (b)  teg  unit. 
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Table  1 

Material  properties  used  in  this  paper. 


Parameters 

Semiconductor  (bismuth-telluride) 

Connector  (Cu) 

X  (W  I(-1  ITT1 

)  /<p  =  /<n 

0.000029T2-0.019593T+  4.809677 

350 

p  (Q  m) 

Pp  —  Pn 

10_6(0.043542T-2.754139) 

1.695  x  10"9 

a  (V  K_1) 

=  ~0Ln 

10"6( -0.002025 T2  +  1 .4234487-44.953611) 

— 

Source 

[27] 

[10] 

N 


7 


sr- 


Fig.  2.  The  effect  of  single  parameter  on  TEG  performance:  (a)  N;  (b)  Hpn;  (c)  y. 


power  servers  as  the  objective  function,  the  maximum  output  po¬ 
wer  can  be  reached  by  looking  for  optimal  parameter  combination 
of  N,  Hpn,  and  y ,  while  the  conversion  efficiency  may  be  reduced 
significantly.  Similarly,  the  maximum  conversion  efficiency  is 
reached  at  the  cost  of  the  output  power.  Therefore,  two  objectives 
namely  maximizing  the  output  power  and  maximizing  the  con¬ 
version  efficiency  of  the  TEG  module  are  considered  simulta¬ 
neously  for  multi-objective  optimization.  In  present  optimization, 
the  multi-objective  function  is  defined  as  follows: 


j 


+  (i  -/) 


V 

Vopt 


(1) 


where, /is  a  weight  factor  for  the  first  objective  function  P,  Popt  is  the 
maximum  value  of  the  objective  functions  P  and  ^0pt  is  the 
maximum  value  of  the  objective  function  77,  when  these  objectives 
are  considered  independently. 


2.3.  Constraint  condition 


Before  the  optimization,  it  is  necessary  to  give  the  reasonable 
range  of  search  variables  based  on  the  practical  requirements. 
Although  single  parameter  analysis  revealed  that  a  small  Hpn  and  a 
large  7  can  improve  P  of  the  TEG  module,  due  to  the  limitation  of 
the  actual  processing,  Hpn  cannot  be  arbitrarily  small,  and  7  cannot 
be  equal  to  1.  In  addition,  larger  N  means  smaller  TEG  unit,  which 
will  bring  processing  difficulties,  therefore,  the  constraint  condition 
of  the  optimization  are  chosen  as  N  <  200,  Hpn  >  0.1  mm,  and 
7  <  0.95. 

2.4.  Optimization  procedure 

Fig.  3  shows  the  flow  chart  of  optimization  procedure.  For 
convenience,  the  search  variables  Ci,  C2  and  C3  are  used  to  denote 
the  semiconductor  pair  number  N,  the  leg  length  Hpn.  The  proce¬ 
dure  for  applying  the  simplified  conjugated-gradient  algorithm  is 
described  briefly  as  follows.  Step  1 :  Make  an  initial  guess  for  the 
search  variables  C\  ~  C3  and  the  values  of  the  search  step  sizes 
0i  ~  fe  with  initialization  accomplished,  the  run  itself  can  begin. 
Step  2:  Use  the  direct  problem  solver  to  predict  the  TEG  module 
performance,  and  calculate  the  objective  function  J.  Step  3:  When 
the  objective  function  reaches  a  maximum,  the  solution  process  is 
terminated,  otherwise,  proceed  to  step  4.  Step  4:  Determine  the 
gradient  functions,  dJ/dCk  (/<  =  1,  2,  3).  Step  5:  Calculate  the 
conjugate-gradient  coefficients^,  and  the  search  directions, 
for  each  of  the  search  variables.  Step  6:  Assign  a  fixed  value  to  the 
step  sizes  fa  (/<  =  1,  2,  3)  for  all  the  search  variables,  and  leave  it 
unchanged  during  the  iteration.  Step7:  Update  the  search  variables 
and  go  back  to  step  2. 

3.  Validation  of  the  optimization  approach 

It  is  well  known  that,  similar  with  other  gradient  algorithms,  the 
simplified  conjugated-gradient  method  may  not  find  the  global 
minimum  of  the  objective  function,  it  can  get  trapped  in  local 
minimum  solutions.  The  solution  that  the  simplified  conjugated- 
gradient  method  arrives  at  depends  on  the  configuration  of  the 
initial  design.  Thus,  several  restarts  may  be  required  to  find  the 
global  solution.  Based  on  the  reason  above,  three  sets  of  initial 
values  of  search  variables,  (N  =  40,  Hpn  =  0.5, 7  =  0.80),  (60, 0.4  mm, 
0.70)  and  (80,  0.6  mm,  0.60),  are  used  in  our  optimization.  Taking  P 
as  the  single-objective  function,  the  multi-parameter  optimization 
is  carried  out  and  the  results  are  shown  in  Fig.  4.  The  output  powers 
are  0.16  W,  0.23  W  and  0.18  W  for  the  three  initial  designs,  then 
they  gradually  increase  as  the  TEG  geometry  is  optimized,  and 
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Fig.  3.  Flow  chart  of  optimization  procedure. 


finally  reach  the  same  maximum  value  of  Pmax  =  0.53  W  at  search 
step  96,  70,  and  64.  The  results  show  that  the  initial  values  of  each 
search  variable  for  the  three  initial  designs  differ  significantly, 
however,  the  same  optimal  design  is  reached,  which  indicates  that 
the  obtained  optimal  design  may  be  a  global  maximum  solution  for 
the  objective  function.  Compared  with  the  initial  designs,  the 
output  power  rises  by  231.25%,  130.43%,  and  194.44%,  respectively. 
Fig.  4b  shows  that  the  conversion  efficiency  for  the  optimal  design 
is  0.72%,  which  is  much  lower  than  that  of  the  three  initial  designs 
(2.46%,  3.00%  and  3.67%).  This  again  confirms  our  analysis  above, 
when  using  a  single  objective  function,  it  is  difficult  to  simulta¬ 
neously  maximize  the  output  power  and  the  conversion  efficiency. 

4.  Results  and  discussion 

A  micro-TEG  module  with  the  base  area  of  only  4x4  mm2  is 
selected  as  the  testing  example  to  display  applicability  of  the  pre¬ 
sent  multi-objective  and  multi-parameter  optimization  approach 


for  TEG  geometry  design.  Due  to  the  dependence  of  P  and  tj  on  the 
applied  current,  it  can  be  expected  that  the  optimal  geometry  of  the 
TEG  module  should  vary  at  different  load  currents.  Therefore,  the 
search  variables  with  initial  values  of  N  =  40,  Hpn  =  0.50  mm, 
7  =  0.80  are  optimized  simultaneously  to  look  for  the  optimal  ge¬ 
ometry  at  AT  =  150  K  with  load  currents  of  J  =  0.10  A,  0.25  A,  0.30  A, 
and  0.50  A.  The  selected  load  currents  represent  the  typical  small, 
medium,  and  large  currents.  It  should  be  pointed  that  the  initial 
design  has  the  maximum  rj  at  /  =  0.25  A  and  the  maximum  P  at 
/  =  0.30  A. 


4.1.  Single-objective  optimization 

To  determine  I^opt  and  ?7opt  in  the  multi-objective  function  (Eq. 
(1)),  two  single-objective  optimizations  with  P  and  p  as  the 
objective  function  are  first  implemented.  Fig.  5  shows  the  evolu¬ 
tions  of  7]  and  P  during  optimization  when  taking  rj  as  the  objective 
function.  For  the  initial  design,  p  is  2.46%,  3.79%,  3.74%,  and  1.61% 
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Fig.  4.  Variations  of  TEG  performance  during  optimization  using  the  three  initial  de¬ 
signs:  (a)  P;  (b)  rj . 


at  I  =  0.10  A,  0.25  A,  0.30  A,  and  I  =  0.50  A.  When  the  TEG  module 
reaches  the  optimal  design,  7]  respectively  increases  to  3.84%, 
3.86%,  3.85%,  and  3.85%.  For  /  =  0.10  A  and  0.50  A,  r/  is  improved 
significantly,  about  56.10%  and  139.13%  higher  than  those  of  the 
initial  design.  However,  for  I  =  0.25  A  and  0.30  A,  p  increases  only 
by  1.85%  and  2.94%.  This  is  mainly  due  to  the  following  two  rea¬ 
sons.  First,  when  thermoelectric  materials  and  temperature  dif¬ 
ference  are  specified,  the  difference  between  maximum 
conversion  efficiencies  for  TEG  modules  with  various  designs  is 
relatively  small.  Second,  the  initial  design  of  the  TEG  module 
reaches  or  approaches  the  maximum  p  at  these  two  currents,  so 
that  the  space  for  the  improvement  of  p  is  small.  Fig.  5b  shows 
that  when  taking  p  as  the  objective  function,  P  of  the  optimal 
design  presents  two  different  change  trends  depending  on  the 
load  current.  At  small  and  medium  currents  of  /  =  0.10  A,  0.25  A 
and  0.30  A,  P  of  the  optimal  design  is  0.13  W,  0.17  W,  and  0.19  W, 
lower  than  0.16  W,  0.28  W,  and  0.30  W  of  the  initial  design. 
However,  at  the  large  current  of  I  =  0.50  A,  P  of  the  optimal  design 
is  0.20  W,  high  than  0.14  W  of  the  initial  design.  This  means  that  it 
is  possible  to  simultaneously  increase  7]  and  P,  and  hence  the 
multi-objective  optimization  is  necessary  and  achievable.  Fig.  6 
shows  evolutions  of  the  three  search  variables  during  optimiza¬ 
tion.  The  optimal  design  depends  on  the  load  current.  The  optimal 
design  has  N  =  48,  Hpn  =  0.74  mm,  y  =  0.55  at  /  =  0.10  A;  N  =  25, 
Hpn  =  0.83  mm,  y  =  0.77  at  I  =  0.25  A;  N  =  23,  Hpn  =  0.73  mm, 


y  =  0.76  at  I  =  0.30  A;  and  N  =  15,  Hpn  =  0.75  mm,  y  =  0.84  at 
/  =  0.5  A. 

Fig.  7  shows  evolutions  of  P  and  p  during  optimization  when 
taking  P  as  the  objective  function.  The  output  power,  P,  is  signifi¬ 
cantly  improved  as  the  three  search  variables  are  optimized  and 
finally  reaches  the  maximum  value  (Fig.  7a),  while  rj  either  in¬ 
creases  or  decreases  depending  on  the  load  current  (Fig.  7b).  The 
output  power  of  the  initial  design  is  0.16  W,  0.28  W,  0.29  W,  0.14  W 
at  I  =  0.10  A,  0.25  A,  0.30  A,  and  0.50  A,  which  is  elevated  to  0.59  W, 
1.16  W,  1.27  W,  and  1.39  W,  respectively.  At  small  and  medium 
currents,  as  expected,  7]  of  the  optimal  design  is  lower  than  that  of 
the  initial  design,  however,  at  large  current  of  0.5  A,  P  and  rj  are 
improved  simultaneously.  It  worth  noting  that  the  improvement  of 
P  is  more  significant  at  higher  currents,  P  of  the  optimal  design  at 
I  =  0.50  A  exceeds  that  at  I  =  0.3  A,  although  the  initial  design  has 
maximum  P  at  /  =  0.30  A. 

4.2.  Multi- objective  optimization 

The  single-objective  optimization  usually  improves  the  output 
power  at  the  cost  of  the  conversion  efficiency  and  vice  versa, 
therefore,  the  multi-objective  optimization  is  implemented  with 
aim  to  enhance  P  and  rj  simultaneously.  Note  that  the  value  of 
weight  factor /in  Eq.  (1)  may  vary  between  0  and  1,  and  any  value 
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gives  the  relative  importance  to  the  single-objective  function.  The 
value  of  weight  factor  can  be  decided  by  the  designer,  and  a  value  of 
0.5  is  adopted  here. 

Fig.  8  shows  the  evolutions  of  the  multi-objective  function,  J, 
during  optimization.  The  initial  value  of  J  is  0.46, 0.61, 0.60,  and  0.26 
at  I  =  0.10  A,  0.25  A,  0.30  A,  and  0.50  A.  As  the  geometry  of  the  TEG 
module  is  optimized  J  increases  and  finally  reaches  0.81, 0.85, 0.89, 


10  20  30  40  50  60  70  80 


search  step 


Fig.  6.  Variations  of  search  variables  during  optimization  of  rj:  (a)  N;  (b)  Hpn;  (c)  y. 


Fig.  8.  Variations  of  the  multi-objective  function  (J). 
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and  0.85,  respectively.  According  to  Eq.  (1),  a  large  J  means  that  P 
and  ri  are  closer  to  their  maximum  value  obtained  in  the  single¬ 
objective  optimization,  however,  J  always  cannot  be  equal  to  1 
unless  the  weight  factor  are  taken  as  0  or  1. 

Fig.  9  shows  evolutions  of  P  and  77  of  the  TEG  module  during 
multi-objective  optimization.  The  output  power  of  the  optimal 
design  is  found  to  be  0.39  W,  1.04  W,  1.27  W,  and  1.16  W,  and  the 
conversion  efficiency  of  the  optimal  is  3.64%,  3.10%,  3.02%,  and 
3.29%.  In  single-objective  optimization,  the  output  power  of  the 
optimal  design  is  0.13  W,  0.17  W,  0.19  W,  and  0.20  W  when  taking 
the  conversion  efficiency  as  the  objective  function,  while  the  con¬ 
version  efficiency  is  1.53%,  2.82%,  3.01%,  and  3.10%  when  taking  the 
output  power  as  the  objective  function.  Comparison  between  these 
data  indicates  that  the  multi-objective  optimization  makes  a 
reasonable  balance  between  P  and  77  to  improve  the  both  simulta¬ 
neously.  Fig.  10  shows  evolutions  of  the  three  search  variables 
during  multi-objective  optimization,  the  optimal  design  has 
N  =  125,  Hpn  =  0.40  mm,  y  =  0.95  at  I  =  0.10  A;  N  =  137, 
Hpn  =  0.13  mm,  y  =  0.95  at  I  =  0.25  A;  N  =  150,  Hpn  =  0.10  mm, 
y  =  0.95  at  /  =  0.30  A;  and  N  =  93,  Hpn  =  0.13  mm,  y  =  0.95  at 
I  =  0.5  A.  It  is  noted  that  the  optimal  designs  for  the  multi-objective 
optimization  require  a  large  N  and  a  small  Hpn,  which  are  signifi¬ 
cantly  different  from  those  for  the  single-objective  optimization 
(Fig.  6). 

Fig.  11  compares  the  I—P  and  1—7]  curves  of  the  initial  and 
optimal  TEG  modules.  The  solid  circles  in  Fig.  11  represent  P  or  77  of 


the  optimal  design,  while  the  hollow  circles  are  their  initial  values. 
It  can  be  seen  that  the  multi-objective  optimization  not  only 
significantly  elevates  P,  but  also  improves  77.  The  interval  of  effective 
load  current,  which  is  defined  as  the  current  interval  when  P  >  0, 
changes  after  optimization.  For  example,  the  maximum  effective 
load  current  for  the  initial  design  is  0.58  A,  it  decreases  to  0.27  A  for 


Fig.  10.  Variations  of  search  variables  during  multi-objective  optimization:  (a)  N;  (b) 
Hp„; (c)  y. 
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the  optimal  design  obtained  at  /  =  0.1  A,  while  increases  to  1.16  A  for 
the  optimal  design  obtained  at  /  =  0.5  A. 

5.  Conclusions 

This  paper  develops  an  optimization  approach  to  design  the 
optimal  structure  of  TEG  module.  In  the  present  optimization 
approach,  the  direct  problem  solver  adopts  a  multi-physics  TEG 
numerical  model  which  has  been  proved  to  be  able  provide  enough 
accurate  prediction  of  the  TEG  performance,  the  inverse  problem 
solver  uses  a  SCG  algorithm.  The  optimization  approach  is  tested  for 
a  bismuth-telluride-based  TEG  module  with  4x4  mm2  base  area. 
The  single-objective  and  multi-objective  optimizations  are 
respectively  implemented  to  look  for  the  optimal  TEG  performance. 
Three  geometric  parameters,  semiconductor  pair  number  N,  leg 
length  Hpn,  and  base  area  ratio  of  semiconductor  columns  to  TEG 
module  y,  are  taken  as  search  variables,  which  are  optimized 
simultaneously. 

When  taking  the  output  power  P  as  the  objective  function,  P  is 
elevated  significantly  by  optimization  of  the  three  geometric  pa¬ 
rameters.  The  output  power  of  the  optimal  design  is  found  to  be 
0.59  W,  1.16  W,  1.27  W  at  /  =  0.10  A,  0.25  A,  0.30  A,  and  0.50  A,  about 
269%,  314%,  338%,  and  893%  higher  than  that  of  the  initial  design. 
However,  the  improvement  of  P  causes  the  remarkable  reduction  in 


1(A) 


1(A) 


Fig.  11.  Comparison  of  TEG  performance  between  the  initial  and  optimal  designs:  (a)  P; 
(b)  V. 


the  conversion  efficiency  77.  Similarly,  when  taking  77  as  the  objec¬ 
tive  function,  the  increase  in  77  for  the  optimal  design  also  accom¬ 
panies  the  significant  reduction  in  P.  Thus,  the  multi-objective 
optimization  is  necessary  to  improve  P  and  77  simultaneously.  The 
present  results  confirm  the  effectiveness  of  the  multi-objective 
optimization.  With  a  weight  factor  of  0.5,  the  optimal  design  has 
P  =  0.39  W  and  77  =  3.64%  at  /  =  0.10  A;  P  =  1.04  W  and  77  =  3.10%  at 
/  =  0.25  A;  P  =  1.27  W  and  77  =  3.02%  at  I  =  0.30  A;  P  =  1.16  W  and 
77  =  3.29%  at  1  =  0.50  A,  which  makes  a  reasonable  balance  between 
P  and  77  to  improve  the  both  simultaneously. 
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